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We report on the anisotropic response, the charge and lattice dynamics of normal and charge- 
ordered phases with horizontal stripes in single crystals of the organic conductor a-(BEDT-TTF)2l3 
determined by dc resistivity, dielectric and optical spectroscopy. An overdamped Drude response 
and a small conductivity anisotropy observed in optics is consistent with a weakly temperature 
dependent dc conductivity and anisotropy at high temperatures. The splitting of the molecular 
vibrations U27 (B u ) evidences the abrupt onset of static charge order below Too = 136 K. The drop of 
optical conductivity measured within the ab plane of the crystal is characterized by an isotropic gap 
that opens of approximately 75 meV with several phonons becoming pronounced below. Conversely, 
the dc conductivity anisotropy rises steeply, attaining at 50 K a value 25 times larger than at 
high temperatures. The dielectric response within this plane reveals two broad relaxation modes of 
strength Ae L D « 5000 and Ae S D « 400, centered at 1 kHz < v w < 100 MHz and v S u ~ 1 MHz. The 
anisotropy of the large-mode (LD) mean relaxation time closely follows the temperature behavior of 
the respective dc conductivity ratio. We argue that this phason-like excitation is best described as 
a long-wavelength excitation of a 2fcp bond-charge density wave expected theoretically for layered 
quarter-filled electronic systems with horizontal stripes. Conversely, based on the theoretically 
expected ferroelectric-like nature of the charge-ordered phase, we associate the small-mode (SD) 
relaxation with the motion of domain-wall pairs, created at the interface between two types of 
domains, along the a and b axes. We also consider other possible theoretical interpretations and 
discuss their limitations. 

PACS numbers: 71.45.-d, 77.22. Gm, 71.30.+h, 78.30.-j 



I. INTRODUCTION 

The competition between the tendency of electrons 
to delocalize and the presence of electron-electron in- 
teraction, electron-phonon interaction - often combined 
with spin-phonon interaction - is the origin of the ex- 
tremely rich phase diagrams in the condensed matter 
systems with reduced dimension, 1 ' 2 Among the most in- 
triguing phenomena found in these systems are broken- 
symmetry phases like charge- and spin-density waves 
(CDW, SDW), charge orders (CO), antiferromagnetic, 
spin-Peierls phases and superconductivity. These phases 
show a large variety of nonlinear properties and complex 
dynamics, including collective excitations. 3-7 . While the 
conventional (2kp) CDW implies periodically modulated 
charge density and is caused by the electron-phonon in- 
teraction, the 4fcp charge modulation of Wigner- crystal 
type, often called charge ordering, originates from strong 
on-site Coulomb repulsion U. Thus, the CO is most com- 
monly thought of as the alternating localized charge of 
different valencies arranged in a crystal lattice. In par- 
ticular, a wealth of CO phenomena with charge dispro- 
portionation is found in organic conductors characterized 
by different anisotropic networks. These include quasi- 
one-dimensional (TMTTF) 2 A: (Ref.i), (DI-DCNQI) 2 Ag 
(Ref. M) and quasi-two-dimensional conductors based on 
the BEDT-TTF [bis(ethylenedithio)tetrathiafulvalene] 
molecule: such as 6>-(BEDT-TTF) 2 RbZn(SCN) 4 and a- 



(BEDT-TTF) 2 I 3 ^ It is notable that these systems have 
a quarter-filled conduction band and due to that not even 
the large values of on-site repulsion U with respect to the 
electron hopping t are sufficient to transform the ground 
state from metallic to insulating. Rather, the inter-site 
Coulomb interaction V is required to stablize a Wigner- 
crystal type phased Although CO features seem to point 
directly toward full localization of charges and their al- 
ternating arrangements, some theoretical studies dispute 
this conventional view indicating that a delocalized CDW 
picture might also be relevant to some CO, e.g., like the 
one in #-(BEDT-TTF) 2 AT materials^ 

a-(BEDT-TTF) 2 I 3 , the first organic material to show 
highly conductive properties in two dimensions^ is 
also one of the most prominent examples of charge or- 
der among 2D organic conductors. This system dis- 
plays a rich temperature-pressure phase diagram with 
a number of intriguing quantum phenomena ranging 
from superconductivity^, charge order^ persistent 
photoconduction, 15 photoinduced phase transition, 16 
nonlinear ultrafast optical response, 17 and zero-gap 
semiconductivityi^ characterized by Dirac cones and 
massless Dirac fermions i 19 : 20 In addition, recently a first 
successful observation was reported of correlated elec- 
tron motion in the light-matter interaction leading to the 
photo-induced insulator-to- metal transition. 21 

The triclinic crystal structure is an alternation of in- 
sulating anion (Fg") layers and conducting layers of donor 
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Figure 1: (a) Schematic representation of donor layer in a- 
(BEDT-TTF) 2 I 3 . Molecular sites belonging to the stack I 
and stack II are denoted as A, A' and B, C, respectively, 
(b) Stripe arrangement in the ab conducting donor layer of 
a-(BEDT-TTF) 2 l3 in the charge- ordered state. Dark- and 
light-gray ovals denote charge-rich and charge-poor molecules, 
respectively. The stripes extend along the ^-direction. 

molecules (BEDT-TTF 5 + on average). The BEDT- 
TTF molecules form a herring-bone structure and are 
organized in a triangular lattice with two types of stacks. 
At room temperature stack I is weakly dimcrized and 
composed of crystallographically equivalent molecules A 
and A', while the stack II is a uniform chain composed of 
B and C molecules [see Fig.QJa)]. Thus the unit cell con- 
tains four BEDT-TTF molecules. At high temperatures 
the system is a semimetal with small electron and hole 
pockets in the Fermi surface! 22 ' 23 A slight but noticable 
charge disproportionation is observed already at room 
temperature, indicating that CO gradually develops as 
the temperature is reduced towards TcO' 24 ' 26 - As demon- 
strated by nuclear magnetic resonance (NMR)— and syn- 
chrotron x-ray diffraction measurements^ charge order 
at long length scales develops fully below the metal-to- 
insulator phase transition Tqo = 136 K. At Tqo the 
conductivity drops by several orders of magnitude and 
a temperature-dependent gap opens in charge and spin 
sector which indicates insulating and diamagnetic na- 
ture of the ground state, x-ray diffraction measure- 
ments indicate subtle structural changes at Ico There 
are no translations of molecules, only a shift in dihe- 
dral angles (angles between molecules of two neighbor- 
ing stacks) is observed. This results in breaking of in- 
version symmetry between A and A' sites, with space 
group changing from PI to the PI, and also allows for 
crystal twinning in the low-temperature acentric struc- 
ture. The variations in dihedral angles cause an appre- 
ciable 2D modulation of overlap integrals between the 
BEDT-TTF sites<2£ Finally, molecular deformations are 
observed which are at the origin of charge disproportion- 
ation. Estimated charge values p of the molecules A, A', 
B and C are 0.82(9)e, 0.29(9)e, 0.73(9)e and 0.26(9)e, re- 
spectively. The exact site assignment is not settled yet, 
since these values differ slightly from those found in the 
NMR, vibrational infrared and Raman spectroscopy, and 
anomalous x-ray diffraction measurements j2£t— Never- 
theless, all these experiments consistently indicate that 
the charge order comprises "horizontal" charge stripes 
of charge-poor (CP) sites, the A' and C molecules, and 



charge-rich (CR) sites (A and B molecules) along the b 
crystallographic axis, as depicted in Fig. [ljb). Contrary 
to the conventional view of the CO as an alternation of 
localized charges, experimental x-ray data indicate that 
a delocalized CDW-like picture might be a more appro- 
priate description of the CO stabilized in this system. 

There have been some previous attempts to find evi- 
dence of collective response in a-(BEDT-TTF) 2 I 3 . More 
than fifteen years ago, preliminary measurements 3 -^ in- 
dicated the existence of a broad relaxation in radio- 
frequency range with a large dielectric constant of the 
order of 10 5 as well as sample-dependent nonlinearities. 
However, no conclusive interpretation could be given at 
that time. These experiments need to be revisited and 
more detailed ones with electric fields along different 
crystallographic directions have to be added in order to 
shed more light on the electrodynamics of the charge- 
ordered state. Recently anisotropic voltage oscillations 
have been observed that are associated with nonlinear 
conductivity^ Such a phenomenon closely resembles the 
collective sliding motion of density waves in quasi-one- 
dimensional conductors^ which makes it yet another re- 
sult to raise the question of similarities and differences 
between conventional density waves and CO found in 2D 
systems. At the time of writing of this paper, Yue et 
al££ presented a new detailed infrared and Raman study 
in the metallic and CO phase of a-(BEDT-TTF) 2 I 3 . 

In an attempt to clarify several issues mentioned 
above we have undertaken an investigation of electrody- 
namic response in the normal phase and ground state 
of the layered organic conductor a-(BEDT-TTF) 2 I 3 . 
We have performed comprehensive optical investiga- 
tions in all three crystallographic directions, as well 
as dc and ac conductivity-anisotropy measurements on 
carefully oriented single crystals of a-(BEDT-TTF) 2 I 3 . 
A brief report of this investigation has previously 
been published? 3 ^ Our results demonstrate a complex 
and anisotropic dispersion in the charge- ordered state 
in contrast to an almost isotropic and temperature- 
independent charge response at high temperatures. The 
abrupt onset of static charge order below Tqo = 
136 K is signaled by a splitting of molecular vibrations 
^27 (B u ) and a dramatic drop of the optical conductiv- 
ity. Its anisotropy does not change much in contrast to 
anisotropy of dc conductivity which rises steeply with 
decreasing temperature. Charge redistribution detected 
in the CO phase corresponds nicely to the one esti- 
mated from the x-ray data. 26 Similar to the Peierls CDW 
state, we observe long- wavelength charge excitations with 
an anisotropic phason-like dispersion, which surface as 
broad screened relaxation modes along the a- and 6-axes 
of the BEDT-TTF planes. In addition, we detect short- 
wavelength charge excitations in the form of domain-wall 
pairs, created due to inversion symmetry breaking, which 
are less mobile and induce a much weaker polarization, 
again along both crystallographic axes. Possible theoreti- 
cal interpretations are discussed and arguments are given 
that the nature of the horizontally-striped CO phase is a 
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cooperative bond-charge density wave rather than a fully 
localized Wigner crystal. 



II. SAMPLE CHARACTERIZATION AND 
EXPERIMENTAL METHODS 

DC resistivity was measured between room tempera- 
ture and 40 K by standard four-contact technique. In 
the frequency range 0.01 Hz-10 MHz the spectra of 
the complex dielectric function were obtained from the 
two-contact complex conductance measured by two set- 
ups. The low- frequency, high-impedance setup covers 
the 0.01 Hz - 3 kHz range, ac voltage signal is applied 
to the sample, the current response is transformed to 
voltage by a Stanford Research Systems SR570 current 
preamplifier and detected using a dual-channel digital 
lock-in Stanford Research Systems SR830. At higher fre- 
quencies (40 Hz - 6 MHz) an Agilent 4294A impedance 
analyzer with virtual ground method was used. Even 
though the impedance analyzer reaches frequencies up 
to 110 MHz, we are limited to approximately 6 MHz 
by cable length. The applied ac signal levels, typically 
50 mV (1 V/cm), were well within the linear response 
regime (upper bound verified to be at least 6 V/cm). 
The results obtained by both methods agree in the over- 
lapping frequency range. In dielectric measurements a 
background contribution of stray capacitances is always 
present due to cabling and sample holder construction. 
In order to account for and remove these influences, we 
have routinely subtracted the open-circuit admittance 
from all measured sample admittances. The background 
capacitance of our setup amounts to 350 fF at all mea- 
sured temperatures. At frequencies 10-10000 cm -1 the 
complex dielectric function was obtained by a Kramers- 
Kronig analysis of the temperature-dependent infrared 
reflectivity measured by Fourier transform spectroscopy 
applying common procedure^ 

The samples under investigation are flat, planar high- 
quality single crystals. As a rule, the pronounced sample 
surface is in the ab plane of the crystal structure. The 
c-axis of the crystal corresponds to the direction perpen- 
dicular to sample plane. The reflectivity measurements 
were performed on as-grown surfaces; for our c-axis in- 
vestigation we employed an infrared microscope. 

Contacts for dc and ac conductivity measurements 
were made by applying carbon paint directly to sample 
surface. On our first studied single crystal contacts have 
been prepared parallel to sample edges without prior ori- 
enting the sample. After dc and ac conductivity mea- 
surements, the orientation of contacts was determined 
to be approximately along [110] (diagonal) direction by 
taking the x-ray back-reflection Laue photograph with 
the beam direction perpendicular to the largest facets of 
the sample. OricntExpress3.3 software^ was applied to 
simulate the recorded flat back reflection patterns. Addi- 
tional experiments were conducted on needle-like samples 
cut from one single crystal that was oriented by recording 



the mid-infrared spectra. The needles were cut along in- 
plane axes which ensured that the orientation of electric 
field is along one of the axes. The dc and ac conductiv- 
ity on those samples was measured parallel to the a-axis 
(perpendicular to stripes; along this direction CR and 
CP sites alternate) and parallel to the 6-axis (along the 
stripes). 

With scrutiny we ruled out any influence of extrinsic 
effects in the dielectric spectroscopy measurements, es- 
pecially those due to contact resistance and surface layer 
capacitance. Of the two relaxation modes detected in the 
dielectric response, the smaller one has a temperature- 
independent relaxation time, thus it can be identified 
as an intrinsic property of the sample. On the other 
hand, the larger dielectric mode bears features which 
might indicate a Maxwell- Wagner-like relaxation at the 
contacts^ and requires a thorough verification of contact 
quality. With this in mind, we performed dc resistance 
measurements in the standard four- and two-contact con- 
figurations. Taking into account the difference in geom- 
etry between contacts in these two configurations, the 
four-point resistance can be scaled and subtracted from 
the two-point resistance in order to estimate contact re- 
sistance (i? c ) vs. sample bulk resistance (i? s )- In the 
insulating phase our samples feature a gradually rising 
R c / R s ratio, from 0.1 just below the transition up to the 
order of 10 at temperatures below 50 K. However, in- 
dependent of the changes in R c /Rs ratio, the strength 
of dielectric response remained approximately constant 
and finite throughout the insulating regime proving that 
the observed behavior of measured capacitance is gov- 
erned predominantly by sample bulk response. Only the 
dielectric spectra in the immediate vicinity of Tqo were 
under influence of instrument artefact and therefore dis- 
regarded. Namely, above 135 K the sample conductivity 
becomes greater than the imaginary part of admittance 
2irfC by a factor of 1000 or more in the whole range of 
frequency /, an instrument limit beyond which no mean- 
ingful measurements of capacitance C can be made. 



III. RESULTS AND ANALYSIS 
A. Optics 

Despite the numerous reports on the optical properties 
of a-(BEDT-TTF) 2 I 3) ^'^ 1 .'2&21^ a few aspects of the 
development right at or below the charge-order transition 
are worthwhile to be reconsidered in the present context. 
Here we concentrate on two issues. The first is the con- 
ductivity and reflectivity in the highly-conducting (ab) 
plane at the metal-to-insulator phase transition which 
yield information on the anisotropy and the energy gap. 
The second issue concerns the redistribution of charge 
on the molecular sites that can be monitored via the vi- 
brational features and their evolution on cooling. Vi- 
brational spectroscopy with light polarized E || c is most 
sensitive to characterize charge orderj 28 ' 44 In addition, we 
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Figure 2: (Color online) Optical properties of a-(BEDT- 
TTF)2l3 for different temperatures as indicated. The upper 
panels (a) and (c) show the reflectivity, the lower panels (b) 
and (d) the corresponding conductivity. Note the different 
vertical scales. On the left side (a) and (b) measurements are 
shown with the electric field polarized parallel to the stacks 
(E || a); the panels on the right side display the data for the 
polarization perpendicular to stacks (E || 6). 



present and discuss vibrational features of the BEDT- 
TTF molecule in metallic and the CO state seen with 
light polarized E || a and b. 



1. Electronic contributions 

The experimentally accessible frequency range extends 
from 10 to 5000 crn^ 1 and covers the bands formed by 
the overlapping orbitals of neighboring molecules. Fig. 
[5] shows the optical properties for the two polarizations 
E || a and E || b in the highly conducting plane at dif- 
ferent temperatures above and below Tqo- The opti- 
cal spectra are dominated by a broad band in the mid- 
infrared in both directions that is different in strength by 
about a factor of 2. A shoulder in R(uj), which shows up 
as a pronounced dip in the conductivity spectra around 
1450 cm -1 , is due to the strong electron-molecular vibra- 
tional (emv) coupling of the v§ (A g ) mode^ 

a. Overall behavior and anisotropy. Although the 
reflectivity exhibits a metallic response at ambient tem- 
peratures leading to a finite conductivity, no Drude-like 
response of the quasi-free carriers can be separated from 
the wing of the mid- infrared band [Fig.^b), (d) and Fig. 
[3]. Hence the weakly temperature-dependent conductiv- 
ity above the metal-insulator transition is described by an 
overdamped Drude response, i.e., a small spectral weight 
compared to the large scattering rate. 

Upon decreasing the temperature from 300 K down 
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Figure 3: (Color online) The far-infrared conductivity for (a) 
E || a and (b) E || 6 for different temperatures above and 
below the charge order transition at Too = 136 K. 
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Figure 4: (Color online) Temperature dependence of the opti- 
cal conductivity of q-(BEDT-TTF) 2 I 3 for T < T C o measured 
for the polarization (a) E || a and (b) E || b. In order to il- 
lustrate the development of the charge-order gap, the phonon 
lines have been subtracted to some extent. Dashed line shows 
the linear extrapolation which gives the optical gap value of 
about 600 cm -1 . Arrows denote the anisotropic dc transport 
gap. 



to Tco, the reflectivity slightly increases due to reduced 
phonon scattering. In the CO state, the far-infrared re- 
flectivity drops dramatically and the corresponding op- 
tical conductivity decreases as the energy gap opens in 
the density of states (Fig. [3J . The spectral weight shifts 
to the mid-infrared range where it piles up in a band 
with maxima around 1500 cm -1 for E || a (parallel to 
the stacks) and 2000 cm -1 for E || b (perpendicular to 
the stacks) 4£ When screening by the conducting charge 
carriers is reduced at T < Tco, the Fano-shaped antires- 
onances in the conductivity due the emv coupled molec- 
ular vibrations become even more pronounced and split 
the mid-infrared peaks. 

b. Optical gap. At the metal-insulator phase transi- 
tion, we see an abrupt opening of an optical gap, with the 
conductivity in the overdamped Drude region dropping 
down to very low values as shown in Fig. [3J In order to 
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Figure 5: (Color online) Development of the spectral weight 
SW(oj c ) as a function of cut-off frequency ui c , calculated by 
SW(aj c ) = 8 Jg" ai(co)dLo — LOp = Anne 2 /m for both direc- 
tions of a-(BEDT-TTF) 2 l 3 . 



illustrate the low-temperature electronic behavior more 
clearly, we have fitted the vibrational features by Lorentz 
and Fano curves, and subtracted them from the measured 
spectra. The results for both polarizations are plotted in 
Fig. @] for different temperatures. The drop of cr(w) be- 
low 1000 cm -1 can be extrapolated linearly to obtain a 
gap value 2A « 600 cm -1 for T — > 0, corresponding to 
approximately 75 meV; it is basically identical for both 
polarizations^ It is worth of noting that while the con- 
ductivity for E || a is indeed close to zero at frequencies 
below 600 cm -1 , the conductivity for E || b remains finite 
down to about 400 cm -1 . Thus taking only the range up 
to 800 cm -1 into account, we can extract gap values of 
600 and 400 cm -1 from the linear extrapolation, which 
corresponds rather well to the one extracted from the dc 
conductivity measurements (see Fig. [5]). 

For T < Too spectral weight still moves from the gap 
region to higher frequencies as T is reduced: it piles up 
around 1000 cm -1 and higher (Fig. [5]). Interestingly, not 
only the region of the gap changes, but spectral weight in 
the entire range shifts to higher frequencies. The maxi- 
mum of the mid-infrared band moves up slightly which in 
part can be described to thermal contraction, but mainly 
to the redistribution of spectral weight. 



2. Vibrational features 

a. Charge disproportionation. In order to charac- 
terize the redistribution of charge values on molecular 
sites associated with the charge ordering we followed the 
behavior of the infrared active i>27 (B u ) charge-sensitive 
mode by measuring perpendicular to the conducting 
plane. This vibration is the out-of-phase contraction of 
the C=C double bonds in the BEDT-TTF rings which 
leads to a dipole-moment change parallel to the long axis 
of BEDT-TTF molecule. The frequency of this mode is 
sensitive to the charge population of the molecule and 
known to split upon passing through a charge-ordering 
phase transition ! 27 ' 44 In the metallic state (T > T"co) 
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Figure 6: (Color online) Temperature dependence of the in- 
tramolecular vibrations of the BEDT-TTF molecule mea- 
sured for the perpendicular direction E || c. The curves for 
different temperatures are shifted by 10 (ficm) -1 for clar- 
ity reasons. The ^27 mode becomes very strong right at the 
charge order transition (T = 295, 210, 160, 136.0, 135.9, 
135.8, 135.7, 135.6, 135.5, 135.0, 80, and 10 K). 
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Figure 7: (Color online) Temperature dependence of 
the intramolecular vibrations ^14 (A 9 ) of the BEDT-TTF 
molecule. The curves for different temperatures are shifted 
by 10 (Slcm) -1 for clarity reasons. For the polarizations (a) 
E || a and (b) E || b the ^14 grows and is splits in three distinct 
peaks (T = 120, 90, 60, 50, 40, 30, and 17 K). 



we observe a wide single band at about 1445 cm ; the 
frequency corresponds to an average charge of +0.5e 
per molecule. The charge diproportionation happens 
abruptly at Tco — 136 K: in the CO insulating state the 
mode splits in two pairs of bands at 1415 and 1428 cm -1 , 
and at 1500 and 1505 cm -1 , as demonstrated in the wa- 
terfall plot of Fig. [TJ] The lower-frequency bands cor- 
respond to approximately +0.8 and +0.85e charge on 
the molecule, the upper-frequency modes to +0.2 and 
+0.15e. This charge redistribution remains constant on 
further cooling and is in agreement with the charges es- 
timated by x-ray for the four different sites in the unit 
cell.— The comprehensive infrared and Raman experi- 
ments of Yue et al££ confirm our findings. 

Coming back to the vibrational features; also for E 



6 



polarized parallel a and b, we observe changes in shape 
of some modes at the metal-insulator transition. Due 
to screening, no molecular vibrations can be seen above 
the metal-to- insulator phase transition. Below Too, the 
modes detected in the in-plane spectra are the features 
of totally symmetric vibrations of BEDT-TTF molecule 
emv coupled with electronic charge-transfer transitions 
which have been observed and assigned previously 39 
down to 300 cm -1 . As a result of interaction with 
electronic transition they have a Fano-shape^ an anti- 
resonance at frequencies where they coincide with the 
electronic excitations and an asymmetric peaks shape 
when the electronic feature is separated in frequen- 
cies. Thus, for example the Agfa) feature and several 
symmetric and asymmetric CH3 vibrations^ 9 - at about 
1400 cm -1 not only shows a blue shift together with 
a charge-transfer band in the mid-infrared, but also 
changes shape to become a narrow and slightly asymmet- 
ric band (see Fig. [2]). While the lower- frequency modes 
are only weakly seen for T > Tco > hi the insulating phase 
the spectra we observe (cf. Fig. [3]) all of the A g vibra- 
tions predicted by Meneghetti et al.£^ for instance, the 
^15 (A g ) mode at 260 cm -1 (associated with the defor- 
mation of the outer EDT rings), the viq(A 9 ) mode at 
124 cm -1 (associated with the deformation of the inner 
TTF rings). These bands are very intense only in the 
E || a direction and barely seen in the b polarization. 
This is in agreement with the symmetry breaking, i.e., 
both the intrinsic dimerisation along the stacks (a direc- 
tion) and the stripes formed along the b direction in the 
CO phase as depicted in Fig. Hfb). A more detailed dis- 
cussion of the in-plane vibrational features was recently 
given by Yue et al. based on low-temperature transmis- 
sion experiments^ 

Interesting, that only the 410 cm -1 mode plotted in 
Fig. [7] changes on cooling in the insulating state. The 
band is much wider than the other features in this range 
at temperatures right below the metal-insulator transi- 
tion and continuously narrows as T is reduced. Follow- 
ing Meneghetti et al^ we assign it to the 1^4 (A g ) mode 
which mainly involves the deformation of the outer rings. 

Finally, the strong vibrational feature observed around 
1300 cm -1 (not shown) is assigned to the emv coupled 
2/4 (A g ) mode of the BEDT-TTF molecule. It is sharper 
and more pronounced for E 1 1 a although the overall con- 
ductivity in the mid-infrared is about half compared to 
E || b. Below 1000 cm -1 a large number of molecular 
and lattice vibrations peak out as soon as the screening 
by the conduction electrons is lost. 



B. Transport 

These measurements reveal that the small dc resistivity 
anisotropy known to be present at room temperature,! 3 
Pal Pb ~ 2, pertains to the whole metallic regime and is 
approximately constant down to Too As a new result 
it has to be pointed out that below Tco the anisotropy 
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Figure 8: (Color online) Resistivity (upper panel) and log- 
arithmic resistivity derivative (lower panel) vs. inverse tem- 
perature of q-(BEDT-TTF) 2 I 3 for E || a (red line) and E || b 
(blue line). 



of resistivity, hence also of conductivity, changes signifi- 
cantly with lowering temperature as the resistivity along 
the a-axis rises more steeply than along the 6-axis, and at 
50 K reaches p a /pb = 50 (see Fig. [5]). In our samples, de- 
spite temperature-dependent activation, the anisotropic 
transport gap in the CO phase for E||a and E||6 can be 
estimated to about 2A = 80 meV and 40 meV, respec- 
tively. At a first glance the appearance of an anisotropic 
transport gap seems to be at odds with the isotropic op- 
tical gap extracted from our optical measurements (Fig. 
|4|- It is well known that systems with a complex band 
structure such as a-(BEDT-TTF)2l3 may exhibit quite 
different optical and transport gaps: optical measure- 
ments examines direct transitions between the valence 
and conduction band, while dc transport probes transi- 
tions with the smallest energy difference between the two 
bands. 

We have also characterized dc resistivity of a-(BEDT- 
TTF) 2 l3 in the conducting ab plane, for E || [110], i.e., 
at an angle of approximately 45° to the crystallographic 
axes. Metallic behavior of resistivity is present from room 
temperature down to 156 K. A sharp metal-to-insulator 
transition is confirmed^ at Tco = 136.2 K, which is ap- 
parent in the peak in d(lnp)/d(l/T) with full width at 
half-height 25T co = 1-5 K; 28T co /T co = 0.011 (Fig. 
[2]). Below the transition the resistivity curve rises with 
a temperature-dependent activation indicating that a 
temperature-dependent conductivity gap opens of about 
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Figure 9: Logarithmic resistivity derivative (main panel) 
and resistivity (inset) vs. inverse temperature of q-(BEDT- 
TTF)2l3 for E || [110], i.e., in the diagonal direction of the ab 
plane. 



80 meV. No significant hysteresis in dc resistivity in the 
vicinity of Tqo could be found. 



C. Dielectric Response 

Low-frequency dielectric spectroscopy measurements 
were performed at various temperatures in the semicon- 
ducting phase. Representative spectra for E || [110] are 
shown in Fig. [TUJ Most notably, between 35 K and up to 
75 K two dielectric relaxation modes are discerned. The 
complex dielectric spectra e(ui) can be described by the 
sum of two generalized Debye functions 



e(w) - £oo = 



Ae 



LD 



Ae S D 



1 + (zwto.ld) 



1 + (zwt 0i sd) 



(1) 

where £oo is the high-frequency dielectric constant, Ae 
is the dielectric strength, tq the mean relaxation time 
and 1 — a the symmetric broadening of the relaxation 
time distribution function of the large (LD) and small 
(SD) dielectric mode. The broadening parameter 1 — a 
of both modes is typically 0.70 ± 0.05. The temperature 
dependences of dielectric strengths and mean relaxation 
times are shown in Fig. [TT] The dielectric strength of 
both modes does not change significantly with temper- 
ature (A^ld ~ 5000, Aesd ~ 400). At approximately 
75 K the large dielectric mode overlaps the small mode. 
It is not clear whether the small dielectric mode disap- 
pears at this temperature or is merely obscured by the 
large dielectric mode due to its relative size. However, 
above 100 K, when the large dielectric mode shifts to 
sufficiently high frequencies, there are no indications of 
a smaller mode centered in the range 10 5 -10 6 Hz. Ac- 
cordingly, above 75 K fits to only a single Debye func- 
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Figure 10: (Color online) Double logarithmic plot of the fre- 
quency dependence of the real (e') and imaginary (e ) part 
of the dielectric function in a-(BEDT-TTF)2l3 at represen- 
tative temperatures for E || [110] . Below 75 K two dielectric 
relaxation modes are observed - full lines for 47 K show a 
fit to a sum of two generalized Debye functions; dashed lines 
represent contributions of the two modes. Above 75 K only 
one mode is detected, and the full lines represent fits to single 
generalized Debye functions. 



tion are performed that we identify with the continua- 
tion of the large dielectric mode. All parameters of the 
large mode - such as dielectric strength, mean relaxation 
time, symmetric broadening of the relaxation time dis- 
tribution function - can be extracted in full detail until 
it exits our frequency window at approximately 130 K. 
At temperatures up to 135 K (just below Tqo — 136 K) 
we can determine only the dielectric relaxation strength 
by measuring the capacitance at 1 MHz. 

One of the most intriguing results is that the tempera- 
ture behavior of the mean relaxation time differs greatly 
between the two dielectric modes. The large dielectric 
mode follows a thermally activated behavior similar to 
the dc resistivity, whereas the small dielectric mode is 
almost temperature-independent. This unexpected and 
novel behavior in the charge-ordered phase raised the 
possibility of anisotropic dielectric response. With this 
in mind we have performed another set of dc and ac 
spectroscopy measurements on the needle-shaped sam- 
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Figure 11: Dielectric strength (upper panel) and mean relax- 
ation time with dc resistivity (points and line, respectively, 
lower panel) in a-(BEDT-TTF)2l3 as a function of inverse 
temperature, for E || [110]. 
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Figure 12: (Color online) Dielectric strength (upper panel) 
and mean relaxation time (lower panel) in a-(BEDT-TTF)2l3 
as a function of inverse temperature; full and empty symbols 
represent parameters of the large and small dielectric mode, 
respectively, for E along the a- (red circles) and b-axis (blue 
triangles). Compared to Fig. 111! the relatively large error bars 
are due to a somewhat unfavorable sample geometry which 
results in higher resistances and a smaller capacitive response. 



pies oriented along the a- and 6-axis. 

Low- frequency dielectric spectroscopy for both E||a 
and E||6 orientation yields results comparable to E |j 
[110]: a large mode whose mean relaxation time fol- 
lows dc resistivity, and a small, temperature-independent 
mode noticeable at temperatures below T ss 75 K. The 
fit parameters to model ([TJ are displayed in Fig. [12] as 
a function of inverse temperature. There is no promi- 
nent anisotropy or temperature dependence in dielectric 
strength, and the Ae values of both the large and small 
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Figure 13: (Color online) Broad-band conductivity spectra of 
a-(BEDT-TTF) 2 I 3 for (a) E || a and (b) parallel & at a few 
selected temperatures. Vertical arrows show the CO optical 
gap. The dashed lines are guides for the eye. 



dielectric modes correspond to those of the sample mea- 
sured in E diagonal orientation. However, an evolution 
of anisotropy in to,ld is clearly visible. Figure Q3] shows 
that the newly-found anisotropy in to.ld closely follows 
the dc conductivity anisotropy. 



IV. DISCUSSION 

It is instructive to compare the conductivity of a- 
(BEDT-TTF)2l3 as a function of temperature in the 
wide frequency range: from dc limit up to terahertz 
frequencies. Fig. [13] composes the conductivity spec- 
tra of a-(BEDT-TTF) 2 I 3 from dc, dielectric and opti- 
cal measurements for E || a and b at different temper- 
atures. First, we address the high-temperature phase. 
The Drude term in room-temperature spectra of or- 
ganic conductors is known to be commonly very weak, 
if present at allj^ The optical conductivity of a-(BEDT- 
TTF)2l3 is in accord to this observation, showing both 
along the a- and 6-axes an overdamped Drude response 
at all temperatures above Tco- The absence of a well- 
defined Drude peak in the vicinity of the CO transi- 
tion resembles the behavior reported for the CO insu- 
lator 6»-(BEDT-TTF) 2 RbZn(SCN) 4 4£ Note that other 
a-phase BEDT-TTF conductors, such as a-(BEDT- 
TTF) 2 MHg(SCN) 4 ,— do exhibit a zero-energy peak. 

It is also worth of noting that the optical, dc resistiv- 
ity and low-frequency dielectric measurements give mu- 
tually consistent values for the conductivity anisotropy 
0fc I (J a ~ 2 at all temperatures above Tco ■ The electronic 
part of optical spectra can be compared to calculations of 
the extended Hubbard model for a quarter-filled square 
lattice using Lanczos diagonalization. They predict a 
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band with a maximum at approximately 6t in the charge 
ordered insulating stated which yields t a = 0.03 eV 
and tb = 0.04 eV for the respective directions. The 
values are in reasonable agreement with Hiickel calcu- 
lations performed by Mori et al^ and support the ob- 
served anisotropy in transport and optical properties. 

Next we address the charge-order phase: how it devel- 
ops on cooling, the ground state features and excitations 
observed by applied spectroscopic techniques. The vi- 
brational spectra reveal that the static charge dispropor- 
tionation sets in rather suddenly (Fig. [6| at the temper- 
ature of MI transition and is accompanied by the respec- 
tive change in the optical properties of the conducting 
plane. At high temperatures, a wide single band at about 
1445 cm^ 1 is observed whose frequency corresponds to 
an average charge of +0.5e per molecule. According to 
Yue et alM- this band originates from slow fluctuations 
of the charge distribution at each site reflecting the par- 
tial charge ordering at short length scales as detected in 
NMR and x-ray measurements! 10 i 26 Yue et al. estimate 
the site-charge distribution slightly above transition to be 
+0.6, +0.6 and +0.4e, agreeing well with the x-ray data 
by Kakiuchi et al. The long-range charge diproportiona- 
tion happens abruptly at Tqo = 136 K and remains con- 
stant on further cooling. In the CO insulating state the 
mode splits in two pairs of bands (see Fig. [6]). The lower- 
frequency bands correspond to approximately +0.8 and 
+0.85e charge on the molecule, and the upper-frequency 
modes to +0.2 and +0.15e, which is in agreement with 
charge estimation by x-ray for the four different sites in 
the unit cel\M 

Interestingly though, optical gap and some of the fea- 
tures of outer ring BEDT-TTF vibrations show a con- 
tinuous change on cooling in the CO state, indicat- 
ing that some changes (but not a charge redistribu- 
tion) happen with temperature in the insulating state. 
In contrast to the sharp onset of the V27 (B n ) vibra- 
tion that monitors the static charge order and does not 
change below Tco (see Fig. [5]), we see some develop- 
ment of the gap and some of the emv coupled features 
on cooling in the charge-ordered state. In the region 
where the gap has opened (T < Tco), the conductiv- 
ity drops further and reaches zero only at the lowest 
measured temperature (T = 17 K). For instance, at 
T = 120 K a finite conductivity is found all the way down 
to 200 cm -1 and even below, in accord with previous mi- 
crowave measurements ! 31 ' 43 The optical gap is more or 
less isotropic, in contrast to the pronounced anisotropy 
of the dc gap, which is explainable taking into account 
that different transitions are involved in optics and dc. 
Nevertheless, as mentioned in Section III, there is some 
weak indication that for E || b excitations are possible 
to lower frequencies. The increase of the anisotropy at 
lower temperatures, which was observed in the dc limit, is 
not that clear in the optical data possibly due to strong 
phonon features and the low-conductivity base line. A 
similar dc conductivity anisotropy has been observed in 
the CO phase of (TMTTF) 2 AsF 6 ^ 
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Figure 14: Anisotropy of the large-dielectric-mode mean- 
relaxation time (points) closely follows the temperature be- 
havior of dc conductivity anisotropy (line) in a-(BEDT- 
TTF) 2 I 3 . 



In the charge-ordered phase we observe novel processes 
at lower frequencies, including CDW responses. As soon 
as the CO phase is entered, the low-frequency conductiv- 
ity drops strongly leading to a step in the radio-frequency 
range (see Fig. [T3|) . This corresponds to the broad and 
strongly temperature-dependent CDW relaxation mode 
(visible only below 120 K), which can be clearly seen in 
the spectra of imaginary part of dielectric function (Fig. 
HOP . It is followed by a power law dispersion attributed to 
hopping transport (for more details on hopping conduc- 
tion see Ref . |6|) , that leads to relatively high ac conductiv- 
ity in the microwave and far-infrared region, as compared 
to the conductivity in dc limit (see Fig. [13]). In the mi- 
crowave region, the most prominent feature is the contin- 
uous increase of conductivity with rising frequency, while 
the far-infrared and infrared regions are mainly charac- 
terized by the suppression of the Drude weight, below 
either the CO gap and strong phonon features. Such 
behavior of conductivity is similar to the one observed 
for fully doped ladders in the (Sr,Ca)i4Cu2404i cuprates 
in which CDW is established. Conversely, comparable 
dc and optical conductivities were found in BaVS3 sys- 
tems in which CDW is also observed^ In this way, we 
can classify the latter system as the fully ordered, while 
ladders and a-(BEDT-TTF)2i3 show features known for 
disordered systems. 

The observed ac conductivity data demonstrate a com- 
plex and anisotropic dispersion in the charge-ordered 
state. First, similar to the Peierls CDW state, we observe 
broad screened relaxation (large dielectric) modes along 
diagonal and both a- and 6-axis of BEDT-TTF planes. 
These modes can be interpreted as signatures of long- 
wavelength charge excitations possessing an anisotropic 
phason-like dispersion. In a-(BEDT-TTF) 2 I 3 Kakiuchi 
et al. were the first to suggest a 2kp CDW which forms 
along the zig-zag path CABA'C of large overlap inte- 
grals detected in their x-ray diffraction measurements— . 
However, the presence of a 2/cp modulation of overlap 
integrals along the p\- and ^-directions, ACA'BA and 
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ABA'CA (see Fig. [T5|) . hints at an additional complex- 
ity and makes the zig-zag paths a somewhat arbitrary 
choice. A theoretical model for the related quarter- 
filled #-ET 2 A systems may provide a more appropriate 
scenario^ Clay et al. showed that the CO phase with 
horizontal stripe phase is characterized by a 1100 modu- 
lation of site charges along the two independent p\- and 
p 2 -directions parallel to the larger overlap integral, and 
with a 1010 modulation along the 6-axis perpendicular 
to the stripes. In addition, this CO is accompanied by a 
2k-p modulation, or tetramerization, of the overlap inte- 
grals along the p-directions, the strongest overlap integral 
being 1-1 and the weakest 0-0. Further, a bond dimeriza- 
tion is also present along the molecular stacks. In other 
words, such a CO phase corresponds to a specific modula- 
tion of bonds and site charges, i.e., a combined 2k-p bond- 
CDW along the two BEDT-TTF plane p-directions with 
bond dimerization in stacking direction. As mentioned 
above, an analogous albeit more complex tetramerization 
of overlap integrals does develop along the p-directions of 
a-(BEDT-TTF) 2 I 3 . Using x-ray diffraction data, Kaki- 
uchi et al. calculated overlap integrals between neighbor- 
ing molecules based on the tight-binding approximation 
and a molecular orbital calculation with the extended 
Hiickel mcthodi 26 i 60 As shown in Fig. [T5j along the p 2 - 
direction, ABA'CA, the strongest overlap integral is ob- 
tained between the two charge- rich sites A and B, quite 
alike the model bond order for the ^-material. Also, bond 
dimerization along the stacking o-direction and its pat- 
tern in the ab plane of a-(BEDT-TTF) 2 I 3 agree with the 
model. However, in the pi-direction, ACA'BA, the or- 
der is shifted by one bond: the largest overlap integral 
is between the charge-rich A site and the charge-poor 
C site. Additionally, the overlap integrals are not per- 
fectly 2fcp sine-modulated along each p-direction. While 
these deviations of a-(BEDT-TTF) 2 I 3 bond order from 
the 6*-ET 2 A model should be recognized, they are hardly 
surprising. Indeed, a-(BEDT-TTF) 2 I 3 has a lower sym- 
metry than a ^-structure which might induce slight dif- 
ferences in bond patterns. Also, the overlap integrals 
obtained from x-ray diffraction could somewhat depend 
on the employed method of calculation. This leaves the 
main physical result of the model by Clay et al., forma- 
tion of a bond-CDW within the conducting molecular 
planes, fully applicable and relevant to the case of a- 
(BEDT-TTF) 2 I 3 . It is plausible to look for the origin 
of phason-like dielectric relaxation in such a 2fcp bond- 
CDW. In this case the energy scale of barrier heights is 
close to the single-particle activation energy indicating 
that screening by single carriers responsible for the dc 
transport is effective for this relaxation. The fact that 
the temperature behavior of the to,ld anisotropy closely 
follows the dc conductivity anisotropy has important im- 
plications: while the CDW motion is responsible for the 
dielectric response, the single electron/hole motion along 
the two p-directions, possibly zig-zagging between them, 
is responsible for the observed dc charge transport. 

The observation of a Peierls-like broad screened dielec- 
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Figure 15: Schematic representation of a 2&f bond-CDW in 
a-(BEDT-TTF) 2 I 3 . Triple, double, single and dashed lines 
show relative strengths of overlap integrals, from strongest 
to weakest, along pi- and f^-directions^. Also denoted are 
dimerized bonds along the a-axis in the AA' and BC stacks. 



trie relaxation in 2D represents an important experimen- 
tal result which clearly indicates that the charge order in 
a-(BEDT-TTF) 2 I 3 cannot be considered of fully local- 
ized Wigner type as predicted by a number of theoretical 
models . 11 ' 53 Rather, as we have argued above, the bond- 
CDW delocalized picture appears as the most appropri- 
ate one. It is noteworthy that a similar 2D dispersion 
was previously observed in the CDW phase developed in 
the ladder layers of Sri4Cu 2 404i. 54-56 Calculations based 
on the extended Hubbard model and the t-J model pre- 
dicted sinusoidal charge-density modulation either of 2/cp 
or 4&f type for this system^ - — This result was at first 
met with surprise due to the strongly correlated nature 
of the ladders system. Nevertheless, the prediction was 
confirmed by resonant x-ray diffraction measurements 56 
which indicated that the established CO is not of fully 
localized Wigner type, but that it can be described as a 
sinusoidal, delocalized modulation. 

Further, we comment on the small dielectric mode 
whose features are characteristic of short-wavelength 
charge excitations. We suggest that the origin of this 
relaxation might be in the twinned nature of the CO 
phase due to the inversion symmetry breaking, with one 
domain being (A,B)-rich and the other (A',B)-rich. 26 In- 
deed, a ferroelectric-like character to the charge-ordered 
phase is suggested by bond-charge dimerization along the 
a-axis together with optical second-harmonic generation 
and photoinduced CO meltingi 17 i 61 ' 62 Our data can be 
most naturally attributed to the motion of charged kink- 
type defects - solitons or domain walls in the charge order 
texture. Both domain walls and solitons stand for short 
wavelength excitations; however whereas a soliton is usu- 
ally a one-dimensional object, the domain wall is not di- 
mensionally restricted. Charge neutrality constraint of 
the CO in a-(BEDT-TTF) 2 I 3 (a change of stripes equiv- 
alent to strictly replacing unit cells of one twin type with 
another) suggests two types of solitons and/or domain 
walls. The first one is the domain wall in pairs (a soliton- 
antisoliton pair) between CR and CP stripes along the 
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Figure 16: Two different types of domain wall pairs in 
the charge-ordered phase of a-(BEDT-TTF) 2 I 3 . (A,B)- and 
(A',B)-rich unit cells are symbolically represented as H — or 
— h cells which form CO stripes. For simplicity we omit the 
B and C molecules. Gray thick lines stand for charge-rich 
stripes. Thin black lines denote a domain wall pair. 



b- axis, which we get if we impose the constraint along 
the 6-axis [Fig. [T6fa)]. The second type of domain- wall 
pair is given by applying the constraint along the a-axis 
so that the domain walls' interior contains both charge 
signs [Fig. ITST b)]. The motion of such entities induces 
a displacement current and can therefore be considered 
as the microscopic origin of polarization in the CO state. 
Namely, in the presence of an external electric field per- 
pendicular to the horizontal stripes, E || a, coupling to 
the AA' dipole moments of each unit cell breaks the sym- 
metry between the two orientations of the dipole. Due to 
first-neighbor interactions the AA' dipoles can most eas- 
ily be flipped at the domain wall, causing the wall pairs 
to move. Coupling to a field in the E || b configuration, 
parallel to stripe direction and therefore perpendicular 
to the AA' dipoles, seems to be more troublesome. How- 
ever, one only needs to remember that A and A' also 
interact with B and C sites. At the energetically unfa- 
vorable domain wall the B and C molecules effectively 
couple AA' dipoles to perpendicular external fields and 
allow for solitonic motion along the 6-axis as well. 

Theoretically, a domain wall should be independent 
on position, however in a real crystal it gets pinned to 
defect sites' Such a pinning causes the domain wall 
to sit in a local energy minimum. The broad distribu- 
tion of relaxation times might then be ascribed to a dis- 
tribution of activation energies associated with pinning 
sites. A weak ac electric field induces a dielectric re- 
sponse which can be attributed to the activation between 
different metastable states over energies barriers. These 
metastable states correspond to local changes of charge 
distribution across the length scale of domain wall thick- 
ness. Finally, a nearly temperature-independent mean 
relaxation time indicates that resistive dissipation cannot 
be dominant for domain wall pairs and that the dielectric 
relaxation is governed by low energy barriers. 

At the end, we consider alternative microscopic de- 



scriptions of excitations which could give rise to the di- 
electric relaxation in a-(BEDT-TTF) 2 l3. A model pro- 
posed by Hotta et alM- describes single-particle exci- 
tations in an insulating, horizontally-striped charge or- 
der on a triangular lattice using plaquette processes. 
These correspond to a localized charged particle mov- 
ing between charge-rich and charge-poor stripes. How- 
ever, the model does not lend itself easily to workable 
predictions of transport, dielectric and infrared proper- 
ties. Another single-particle description might be found 
in the excitonic model used by Yamaguchi et al. to 
explain dielectric properties below T — 2 K as well 
as nonlinear conductivity of 6»-BEDT-TTF 2 MZn(SCN) 4 
(M = Cs, Rb)i 65 ' 66 This model considers a charge- 
ordered quarter-filled square-lattice system in which ex- 
citations are created by moving a localized hole to a site 
where no hole was originally present, thus creating a 
pair of a localized electron and hole. Due to their elec- 
tric fields being confined to the polarizable plane of ET 
molecules, the attractive Coulomb potential between a 
bound electron and hole is of a logarithmic form: the 
potential is modeled as U — Uo\n(r/a) for r < A, and 
a constant U = Uq\u{\/ a) for r > A, where the spa- 
tial scale a is the distance between neighboring BEDT- 
TTF sites, a « 0.5 nm, and A is the screening length. 
At finite temperatures bound pairs are thermally ex- 
cited. The potential barrier for a pair to unbind under 
an electric field E is given by 2A(E) w Uq In {Uo/(eEa)}, 
E -c Uo/(ea). In the excitonic picture a transport gap 
2A(E = 0) = Uo In (A/a) naturally follows from the cut- 
off of logarithmic potential^ The thermally excited but 
bound pairs are polarized in the presence of an electric 
field and give rise to a temperature-dependent dielectric 
constant 
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where no is the electron-hole density at T — > oo which we 
take to be equal to half the BEDT-TTF density^ The 
above expression can be compared with the measured di- 
electric properties of a-(BEDT-TTF) 2 I 3 . Namely, the 
total dielectric constant [e(w — 0)] is well-approximated 
by the Ae of LD mode (see Figs. [11] and [T^J) . A two- 
parameter (A, Uq) fit to, e.g., our E || a data above 
T = 25 K reproduces adequately the general temperature 
dependence, but gives a rather small A = 2.2 ±0.1 run 
and Uo — 1.6 ± 1.5 meV. The fit value of Uo is in 
stark contrast with the value extracted from transport 
gap, Uo — A/ln(A/a) « 50 meV (further, substitut- 
ing A from the previous expression gives an unsatisfac- 
tory one-parameter fit). A similar discrepancy can be 
seen in E || [110] and E || 6 fits. Also, it is not clear 
whether this model (here shown only in the static limit) 
reproduces the shapes of our experimental dielectric spec- 
tra and their temperature dependence. In the end the 
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excitonic picture, while nicely applicable to 0-BEDT- 
TTF 2 MZn(SCN) 4 below 2 K, does not seem to account 
for general dielectric features of a-(BEDT-TTF) 2 l3 in 
the charge-ordered phase. 

V. SUMMARY 

We investigated electrodynamic properties in the 
single-crystals of the layered organic compound a- 
(BEDT-TTF^Is- In the normal phase, we observe an 
overdamped Drude response and a weak optical con- 
ductivity anisotropy. This is consistent with an almost 
isotropic, weakly temperature-dependent dc conductivity 
inside the conducting layers. Broad intramolecular vibra- 
tions of the BEDT-TTF molecule might be attributed to 
charge ordering fluctuations which form at short length 
scales. We demonstrate the abrupt onset of static charge 
order below Tco = 136 K followed by a dramatic drop 
of the optical conductivity. The charge diproportiona- 
tion remains constant on further cooling. The observed 
charge values are +0.8 and +0.85e on charge-rich sites, 
and +0.2 and +0.15e on charge-poor sites, consistent 
with the charges estimated by x-ray. Below the charge- 
order transition we detect the strong development of in- 
plane dc conductivity and dc gap anisotropy in contrast 
to the anisotropy of optical conductivity which remains 
weak and similar to high temperatures. The optical gap 
is approximately 75 meV. The development of dc con- 
ductivity anisotropy in the conducting layers is accom- 
panied by appearance of two dielectric relaxation modes 
in kHz-MHz frequency range. The large dielectric mode 



features an anisotropic phason-like behavior, whereas the 
small dielectric mode is temperature-independent and its 
properties are reminiscent of a soliton-like behavior. All 
these results make it rather clear that the most consistent 
picture of the horizontal stripe, charge-ordered state in 
layered a-(BEDT-TTF)2l3 is a cooperative bond-charge 
density wave with ferroelectric-like nature, rather than 
a fully localized Wigner- crystal. Additional theoretical 
and experimental work is warranted in order to refine 
the microscopic description of both metallic and insulat- 
ing phase of a-(BEDT-TTF) 2 l3. The emerging broader 
issue of the anisotropic, screened phason-like dispersion, 
its associated sinusoidal charge-density modulation and 
whether they are a general signature of charge order 
in layered strongly correlated systems certainly deserves 
further efforts. 
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